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A magnetization process in two-dimensional ferrimagnet BIPNNBNO is analyzed. 
The compound consists of ferrimagnetic (1,1/2) chains coupled by two sorts of anti- 
ferromagnetic interactions. Whereas a behavior of the magnetization curve in higher 
magnetic fields can be understood within a process for the separate ferrimagnetic 
chain, an appearance of the singlet plateau at lower fields is an example of non-Lieb- 
Mattis type ferrimagnetism. By using the exact diagonalization technique for a finite 
clusters of sizes 4x8 and 4 x 10 we show that the interchain frustration coupling 
plays an essential role in stabilization of the singlet phase. These results are comple- 
mented by an analysis of four cylindrically coupled ferrimagnetic (1,1/2) chains via 
an abelian bosonization technique and an effective theory based on the XXZ spin- 
1/2 Heisenberg model when the interchain interactions are sufficiently weak/strong, 
respectively. 



I. INTRODUCTION 



During the last fifteen years, two-dimensional (2D) quantum spin systems have attracted 
a lot of attention both from theoretical and experimental physicists. A competition between 
conventional classically ordered phases and more exotic quantum ordered phases lies in 
the focus of the investigations. Magnetic systems with a finite correlation length at zero 
temperature and a finite spin gap above the singlet ground state, spin liquids, realize Haldane 
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prediction at the level of two space dimensionsi. 

To date, one can distinguish two main routes in studies of 2D spin gap compounds. A 
formation of spin gap in spin dimer systems, for example SrCu2(B03)2^ and CaV40g^, is 
explained by a modified exchange topology similar to Shastry- Sutherland lattice^. Another 
way to increase quantum fluctuations and stabilize a spin liquid ground state is realized 
in kagome antiferromagnet^. Experimental candidates for 2D kagome antiferromagnets are 
currently available: herbertsmithite^"- and volborthite^ii^. Both these strategies deal with 
antiferromagnetic compounds. In view of this, an observation of a singlet ground state with 
a pronounced spin gap in 2D ferrimagnetic material BIPNNBNO seems exotic^. 

The crystal structure of BIPNNBNO is shown in Fig. [H A magnetic unit of the spin 
system presents organic triradical BIPNNBNO. Each of the molecules includes three s=l/2 
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FIG. 1: Magnetic model of BIPNNBNO crystal. The black (white) circles denote spins S=l 
(s=l/2). 

spins (see Fig. |2]) with intramolecular ferromagnetic Jp and antiferromagnetic Jaf interac- 
tions. The magnitude of jJi?! ~ 300 K is very large, and two spins coupled ferromagnetically 
behave as a S - 1 moiety. Ferrimagnetic chains are stretched along the b-axis. There are two 
kinds of antiferromagnetic interchain interactions along the a-axis. One is between the s-1/2 
spins, which connects the nearest neighboring chains. The other is between S-1 species, 
which connects the next nearest neighboring chains and introduces spin frustration. 
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FIG. 2: Molecular structure of BIPNNBNO and the elementary magnetic cell at Jp ^ |^af|- 

The puzzle is following. It is well known that a low-energy physics of an isolated ferrimag- 
netic (S,s) chain corresponds to a gapless (S-s) ferromagnelr^. Quite predictably one might 
expect an appearance of an ordered state with fluctuations in the form of spin waves near 
the classical state. However, measurements of magnetization shows that an opening of a gap 
by analogy with the Haldane chain is likely scenario. Such a behavior is a manifestation of 
non-Lieb-Mattis type ferrimagnetism^^. Namely, the magnetization measured at 400 mK is 
nearly zero below 4.5 T, increases rapidly above 4.5 T and exhibits a broad 1/3 plateau and 
a narrow 2/3 one at 7-23 T and around 26 T, respectively. Above 29 T, the magnetization 
is completely saturated^^. 

The purpose of the paper is to investigate the magnetization process. The problem is com- 
plicated by a lack of reliable information about intra- and interchain exchange interactions. 
So, before studying of the 2D ferrimagnetic system BIPNNBNO we develop a simple quan- 
tum mechanical approach that models a magnetization process of the ferrimagnetic chain 
(1,1/2). The treatment agrees qualitatively with a predictions of the theory for quantum 
spin chainsi-^ and provides reasonable estimations of the exchange intrachain couplings. In 
addition, it captures a peculiarity of a magnetization process in the prototype 2D material, 
i.e. an appearance of the intermediate 2/3 plateau. Given these estimations we examine the 
magnetization process in BIPNNBNO by analyzing exact diagonalization (ED) calculations 
for a finite clusters of size = 32 and N = 40. A main conclusion to be drawn from these 
calculations that an emergence of the anomalous singlet plateau is a consequence of the 
frustrating interchain interaction. 

Two different mechanisms of formation of the plateau may be likely candidates: a gener- 
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alization of Haldanes conjecture to the weakly coupled ferrimagnetic chains, and a valence- 
bond (dimerized) type ground state in the strong-coupling limit. To determine what of these 
scenarios is relevant we develop low-energy effective theories for the 4-legs spin tube, which 
forms a minimal setup including the interchain couplings. In the regime of weakly coupled 
spin tube legs we apply abelian bosonization technique. The opposite limit of a strong-ring 
interaction is analyzed in terms of an effective Heisenberg XXZ model where the intrachain 
coupling is perturbatively taken into account. Our analytical treatment shows that only the 
first approach confirms an important role of frustration in stabilization of the singlet phase. 

Note that a study of spin tubes is of interest by itself because both of frustration and quan- 
tum fluctuation are strong^^. Our model is directly related with the compound BIPNNBNO, 
but the main results are expected to apply to other frustrated spin tubes as well. Re- 
cently, it has been reported that the experimental candidate for the four-leg spin tube, 
Cu2Cl4-D8C4S02, is availablei^. 

The paper is organized as follows: In Sec. II we consider a magnetization process in the 
ferrimagnetic chain (1, 1/2). In Sec. Ill we discuss results of magnetization process in two- 
dimensional ferrimagnetic system obtained via the exact diagonalization method on a finite 
cluster. In Sec. IV we derive effective low-energy spin-1/2 Hamiltonian. A bosonisation 
study of the spin tube is carried out in Sec. V. A discussion of these results is relegated to 
the Conclusion part. 

II. MAGNETIZATION OF AN ISOLATED FERRIMAGNETIC CHAIN 

The issue that we address below is whether a calculation for an isolated ferrimagnetic 
chain partially reproduces features of the magnetization curve observed in the BIPNNBNO 
crystal. We demonstrate that both the 1/3 plateau and the 2/3 plateau can be recovered 
within a simple quantum mechanical analysis of a magnetization process of an isolated 
quantum (1,1/2) ferrimagnetic chain under an applied magnetic field. This enables to 
estimate the intrachain exchange parameters Jaf and Ji. 
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FIG. 3: Schematic picture of states of the ferrimagnetic chain used in construction of a magneti- 
zation curve. Excited blocks are marked by the gray shadow. 

We start from the values of the critical fields derived from the ED datap^ 
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- E{N + 1) - E{N), 
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where e, m are the energy and the magnetization per an elementary cell of the (1,1/2) 
ferrimagnetic chain, is a number of the elementary cells. 

To estimate the energies E{S) in the right-hand side of Eqs.([T]), where 5* is a total spin 
of the chain, we use the Hamiltonian of the one-dimensional quantum (1,1/2) ferrimagnet 

TV N 

JaF ^ SiiS2i + -^1 ^ S2iSii+i {Si = 1, S2 = 1/2), (2) 
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and construct the required quantum states \SM) with the given quantum numbers of the 
total spin 5* and its z-projection M. 

We suppose that the 1/3 magnetization plateau corresponds to the ground state of the 
(1,1/2) ferrimagnetic chain with S = N/2. The wave function of the polarized state is given 
by 



N 



\N/2,N/2) = l[ 



i=l 



2/22 



(3) 



and presents a direct product of the spin states of the magnetic elementary cells (1,1/2) [see 
Fig. |3]. The corresponding energy eigenvalue equals to 



EiN/2) = -JafN - -JiN. 



(4) 



With an increasing of a magnetic field the ground state ([3]) is destroying and the state 
with S = N/2 + 1 stabilizes. A low-lying excitation may be qualitatively considered as a 
forming of one triplet bond. The trial new wave function is 
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|iV/2 + l,iV/2 + l) = -^J] 
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J^^k'^k- (5) 



k=l 



It is composed from all arrangements of the excited block within the chain taken with equal 
weights ttfc. 

By introducing the state and calculating the matrix element 



EiN/2) + pAF + ^Ji 
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one obtains the relationship for the coefficients 



E{N/2) + ^-Jaf + ^Ji - E{N/2 + 1) 



■^Jl^k,k'±l 



Oik - -^JlOik+l 



(6) 



(7) 



which is tantamount to 



E{N/2 + 1) = E{N/2) + ^Jaf + ^ Ji - \ji 



(8) 



Oik oik J 

This expression includes two independent variational parameters Uk-i/uk and 0:^+1/0;^. The 



minimal value 



E^UN/'l + 1) = E{N/2) + Ijaf - ^Ji 

Z io 



(9) 
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is reached provided ak-i/o-k = oik+i/oik = 1- This yields the critical magnetic field Bi 
destroying the 1/3 plateau 

Bl = - ^Jl. (10) 

To find the critical fields B2 and -B3 of the beginning and the end of the 2/3 plateau, 
respectively, we construct the trial states 
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which are schematically shown in Fig. [31 
By the same manner we obtain 

The saturation field B^„f is determined with an aid of the trial wave functions 
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(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



Given the experimental estimations for the 2D BIPNNBNO system, Bi ~ 31 K, B2 ~ 
B3 ~ 35 K, and -Bsat = 39 K, we obtain from Eqs. fll0lll7p the values of the intrachain 
exchange couplings, Jaf ^ 21 K and Ji ~ 3.5 K. By substituting them into Eq. (|T4|) we get 
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the critical fields of the 2/3 plateau, B2 ^ 32.8 K (24.4 T) and B3 ^ 34.4 (25.6 T). A 
qualitative behavior of the ferrimagnetic chain magnetization curve built from these reference 
points is depicted in Fig. HI We emphasize especially that an emergence of the intermediate 
2/3 plateau is not related with interchain frustration effects. 




FIG. 4: Qualitative magnetization curve of the ferrimagnetic (1, 1/2) chain. 



III. MAGNETIZATION: EXACT DIAGONALIZATION 



J nn I I I I 



FIG. 5: Cluster of = 32 sites used in the exact diagonalization. The fixed exchange couplings 
are Jaf = ^IK, Ji = 3.5K, Jnn = 0.5Ji = 1.75K. 



In order to understand a role of the interchain couplings the magnetization process of 
the BIPNNBNO ferrimagnet was examined by the variant of a numerical diagonalization 
method with conservation of a total cluster spip^"^ 
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The model Hamiltonian is given by 




of the antiferromagnetic interactions, Jaf^ Ji, Jnn and Jnnn, are shown in Fig. [51 We 
perform calculation of the N-step magnetization curve for the N=32 cluster depicted in the 
same Figure. The intrachain parameters, Jaf and Ji, have been estimated in the previous 
Section whereas the interchain ones, Jnn and Jnnn, are assumed to be less than Ji. The open 
boundary conditions are used for the numerical calculations. 

The magnetization process is compared with the results of the model of non-interacting 
(1,1/2) ferrimagnetic chains. A standard way to build magnetization curve at T = 
is to define the lowest energy E{N, M) of the Hamiltonian ([2]) in the subspace where 
Ylf=i i^j + ■^j) = M for a finite system of N elementary (5*, s) blocks. Applying a magnetic 
field B leads to a Zeeman splitting of the energy levels, and therefore level crossing occurs on 
increasing the field. These level crossing correspond to jumps in the magnetization until the 
fully polarized state is reached at a certain value of the magnetic field. The magnetization 
of four independent chains is then derived from 



which gives a step curve. 

An importance of the frustrating coupling is seen from comparison of two magnetization 
curves displayed in Figs. El [71 They correspond to no frustration case and a pronounced 
frustrating coupling, respectively. The magnetization curves exhibit several interesting fea- 
tures. For instance, the magnetization behavior in higher magnetic fields {B > Bi) is well 
reproduced within the model of non-interacting ferrimagnetic chains. Another remarkable 
feature revealed by Figs. [6l [3 that the singlet ground state plateau emerges at non-zero 
frustration interaction whereas the narrow 2/3 plateau appears regardless of the frustration. 
We numerically found that the width A5 of the singlet plateau scales almost linearly with 
a Jnnn valuc (Fig. [H]). To check into the case of the dependence we repeat cacluations on 
a cluster of larger size, = 40, with the same set of parameters that support the finding. 
The observation points out that the zero magnetization plateau has a quantum origin with 
a crucial role of frustration which destroys a long-range order and drives the system into the 
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FIG. 6: Magnetization curve for the 32-site cluster. The exchange coupHngs are taken as Jaf = 
21K, Ji = 3.5K, Jnn = 0.5Ji, and Jnnn = (no frustration). The dotted hne marks a calculation 
via the model of non-interacting (1, 1/2) chains. 




FIG. 7: Magnetization curve for the 32-site cluster. The exchange couplings are taken as Jaf = 
21K, Ji = 3.5K, Jjin = 0.5Ji, Jnnn = 0.075Ji. The dotted line marks a calculation via the model 
of non-interacting (1,1/2) chains. 
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FIG. 8: Value of the singlet plateau as a function of the frustrating coupling Jnnn obtained on 
a cluster of size = 32 (black circles) and = 40 (white circles). 

singlet phase. Below, we address analytically the issue in the regimes of strong and weak 
interchain couplings. 



In low-dimensional Heisenberg systems frustrating couplings can drive transitions to gap- 
full quantum states, where local singlets form a ground state. These quantum gapped phases 
may have long-ranged singlet order (valence bond state), or realize a resonating valence band 
spin liquid. In last ground state is a coherent superposition of all lattice-coverings 

by local singlets^. 

To recognize features of these phases in the ED results we undertake analytical treatments 
of the four-legs spin tube shown in Fig. [9l The new system is infinite along the 6-axis, 
and periodic with the 4-site period along the a-axis. The tube forms a minimal setup 
including the interchain nearest- and next-to-nearest neighbor couplings and contains the 
same number of ferrimagnetic chains parallel to the b-axis as the clusters in the ED study. 
As we demonstrate below, the simplified model elucidates an important role of frustration 
in stabilization of the singlet phase. 



IV. A FORMATION OF THE SINGLET PLATEAU 
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FIG. 9: 4-leg spin tube structure used in the limit of strong ring coupling (up) and in the case 
of weakly interacting chains (below). The black (white) circles mean spin-1 (spin-1/2) sites. The 
gray circles denote renormalized spin-1/2 blocks. 

The singlet phase may arise in the limit of strong ring coupling Jnn, J mm ^ Ji- In this 
case, the problem can be analyzed in terms of Heisenberg XXZ model similar to ladders in a 
magnetic field™. The opposite limit ( Ji ^ Jnn, -^mm) results in a scenario of weakly interact- 
ing chains. Based on a block renormalization procedure the original system is then mapped 
onto the model of a spin tube with four ferromagnetic spin-1/2 legs. We mention that the 
ground state properties of two-leg spin ladders with ferromagnetic intrachain coupling and 
antiferromagnetic interchain couplings have been discussed in Refs.— in an absence of an 
external field. A magnetization process of those spin ladders with an even number of legs 
(2 and 4) has been studied in Ref.— in the regime of weak ferromagnetic coupling along the 
legs and strong antiferromagnetic coupling along the rungs. 

An appearance of the singlet phase in the frustrated spin tube with four weakly coupled 
ferromagnetic spin-1/2 legs can be studied through the bosonization technique which proves 
effectiveness for quasi-one-dimensional spin-one-half systems. To the best of our knowledge, 
the system has never been previously reported, however our further analysis follows closely 
to that of given in Ref.'^, where 4-legs spin tube with antiferromagnetic chains and a specific 
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form of diagonal rung interactions (but with no frustration) has been treated. Note as well 
that spin ladders with ferromagnetic and ferrimagnetic legs are much less studied^^i^I by the 
bosonization approach in comparison with ladder systems with antiferromagnetic legs. The 
main problem arising here is that the formalism is well defined only if there is an easy-plane 
exchange anisotropy. In this regard, we note that measurements of the angular dependence 
of the ESR linewidth for the BIPNNBNO system showed that the largest linewidth was 
observed for the field direction perpendicular to the ab plane^^. Due to the theoretical 
consideration by Oshikawa and Affleck^- a critical regime of XY-anisotropy is expected in 
the compound. In addition we point out that the ED algorithm invoked in the previous 
Section enables to treat clusters of sufficiently large sizes due to use of the rotational SU{2) 
symmetry. The latter is broken by the anisotropy whose role in a singlet gap formation is a 
subject of future ED studies. 



A. Spin tube: weakly interacting rings and a model of a single XXZ chain 

We study the Hamiltonian of the spin tube (see Fig. [9]) 

N N N 4 

n=l n=l n=l i=l 

(20) 

where the Hamiltonian of the separate ring is 

N 

ffring _ j^^ (s„ iS„ 2 + ^n,2^n,3 + ^n,3^n,4: + ^n^i^n,!) + Jnnn ^ ^ (S„,,lS„,^3 + S„,^2Sn,4) • 

n=l 

Here, 5 = 1 and s = 1/2, n is the index of the ring, is the total number of rings, and the 
index i marks the (1,1/2) blocks inside the rings. Periodic boundary conditions along the 
tube direction are imposed. In our model it is suggested that Jnn, Jnnn ^ Ji- 

In the limit Ji = the system decouples into a collection of nonineracting rings. At zero 
magnetic field the singlet and triplet states 

/3 1 
iV^o) = -^|00;00) + -|ll;00), 

|^i) = i=|01;ll) + i=|10;ll) (21) 

have the lowest energies Eq = —2 7,^/9 + 8 Jnnn /9 and Ei = — Jnn/ 9 + 8 Jnnn /9, respectively. 
The states of the ring IS'125'34; SM) are obtained via the common rule of addition of moments, 
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where 5*12 (5*34) is spin of dimer composed of the spins of the 1 and 2 (3 and 4) blocks. The 
singlet and triplet states of the ring that enter into f l2T]) are given in Appendix A. 

Upon increasing the magnetic field a transition between the singlet and triplet states 
occurs at -B = Jnn/9 and the the total magnetization jumps abruptly from zero to M = A^. 

At non-zero ring coupling the sharp transition is broadened and starts from a critical value 
Bq. To find the field we derive the XXZ spin chain Hamiltonian by we using the standard 
approach that is analogous to study of spin- 1/2 ladder with strong rung exchangers. 

The Hamiltonian ( 120|) is splitted into two parts 

N AN 



n=l i=l n=l 

N AN 

n,lS„,+ l,l + S„^2S„+l,2 + S„,3S„+i_3 + S„_4S„+i,4] — {B — Be) ^ ^('S'j^n + ^in)y 
n=l i=l n=l 

where B^. = Ei — Eq. The "Hi lifts the 2^-fold degeneracy of the ground state of the 
Hamiltonian "Ho • The later can be either in the state \iPq) or {ipi)- By using the standard 
many body perturbation theory2^ the effective Hamiltonian can be derived 

N N N 

nj _ jeff \ ^ I QX QX I oy oy \ , yeff \ ^ qz qz tjcS \ ^ qz fOO\ 

'T-eff — 'Jxy I ^n^n+l \ ^ ^ z z '-'n^ \'^^) 

n=l n=l n=l 

where Jfy = -16Ji/27, Jf = -Ji/9 and Bf = Ji/9 + B - B^. 

To get the expression the pseudo-spin S'j = 1/2 operators that act on the states Itpo) and 
\ipi) are introduced 

Sn\'^o)n = \i^l)n, ^+|V^l)„ = 0, 

S-\^0)n = 0, S-\^i)n = mn. (23) 

The starting spin-1 operators and the pseudo-spin operators in the restricted space are 
related by 

1 1 ~ 

= _ _|_ _ 

Si = i-iy-'^Sl, 5r„ = i-^y-'^S-. (24) 
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The corresponding map for the spin-1/2 operators is 

1 1 ~ 

o2 ^ t 

^in 24 12 

The Jordan- Wigner transformation maps the Hamiltonian ( 122|) onto a system of inter- 
acting spinless fermions 

N N N 

n=l n=l n=l 

where t = J^J/2, \/ = Jf and fx = Jf + Bf. 

The lowest critical field Bo corresponds to that value of the chemical potential fi for which 
the band of spinless fermions starts to fill up. This yields the condition /x = —2t and leads 
to the result 



1 16 

Bo = g Jnn + ^Ji. (27) 

The critical value involves no frustration parameter Jnnn that is clearly contrary to the ED 
results. 

A similar analysis can be carried out for the saturation field. The details of the calcula- 
tions are relegated to Appendix B. 



B. Spin tube: a model of weakly interacting ferromagnetic legs and abelian 

bosonization 

To apply the bosonization we should map the initial system, consisting of two sorts of 
spins, spin-1/2 and spin-1, to the spin-1/2 system by using the quantum renormalization 
group (QRG) in real space based on the block renormalization procedure^^. To exploit 
the real-space QRG technique one divide the spin lattice into small blocks, namely, the 
intrachain dimers (1,1/2), and obtains the lowest energy states {|a)} of each isolated block. 
The effect of inter-block interactions is then taken into account by constructing an effective 
Hamiltonian H^g which now acts on a smaller Hilbert space embedded in the original one. 
In this new Hilbert space each of the former blocks is treated as a single site. The effective 

N 

Hamiltonians H^s = HQ is constructed via the projection operator Q = Y\Qi with 

i=l 
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m 

Qi = ^ \<y) (a I of each i-th block where m is the number of low energy states that are kept 
0=1 

and is a number of lattice cells. 

We hold the lowest doublet S — 1/2 to find an effective low-energy Hamiltonian. The 
higher energy 5* — 3/2 states are neglected. One can check that the reduced matrix elements 
are {Si = 1, S2 = 1/2) 

(ll/2;l/2||5i||ll/2;l/2) = 2y|, 

(ll/2;l/2||s2||ll/2;l/2>=--^. 
Therefore the effective spin-1/2 operators of the renormalized chain are 

QlSuQ^ = Qls2^Q, = -^5, {S = 1/2). (28) 

The renormalized Hamiltonian of the intrachain interactions corresponds to the ferromag- 
netic Heisenberg spin-1/2 model with the exchange coupling J = —AJi/9 (Fig. |9]). The 
interchain interactions between the nearest neighbors, spins -1/2, and next to the nearest 
neighbors, spins -1, are renormalized as J± = Jnn/9 and = 16Jnnn/9, respectively. 

Consider a four-legs spin tube consisting of spin-1/2 chains. The Hamiltonian of the 
system is 

4 

"Htube = "Ha + "^12 + "^^3 -|- ?i^4 -|- ?i^4 -|- "Hi^ + 7^24"- (^9) 

A=l 

The spins along the chains are coupled ferromagnetically, the Hamiltonian for the separate 
A-th chain is 

N N 

nj jxy \ ^ / QX QX j_ cU cU \ jz \ ^ qz qz 

rt\ — J WAj'^AJ+I "T 'JAj^Aj+iJ / . '^Aj^AJ+D 

i=l 1=1 

where S^^'J'^ are the spin S=l/2 operators at the jth site, the intraleg coupling is ferromag- 
netic, J > 0. 

The interaction parts are given by 

N N 
"^AA' = '^±,XX' ('^A.j'S'vj + Sf jSyj) + J± XX' '^AJ'S'aj, (30) 

j=l i=l 

N N 

nj' ± _ l' \^ ( qx qx \ CV qy \ \ j' Z \^ qz qZ ( n-,\ 

n-AA' — -'±,AA' 2-^ WXJ-^X'J + '^A.i'^A'J;' + ±,XX' '^Aj'^Aj \'^^) 
j=l i=l 
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and includes the nearest, J± > 0, and the next-to-nearest, > 0, antiferromagnetic inter leg 
couplings. 

The unitary transformation keeping spin commutation relations 

maps the Hamiltonian fl29l) to the Hamiltonian with antiferromagnetic legs. It changes 
jxy _s. —j^y and — )■ J^, and the ferromagnetic isotropic point is A = J^/J^^ = — 1 in 
the Hamiltonian with the antiferromagnetic legs. 

Following the general procedure of transforming a spin model to an effective model of 
continuum field, we convert the spin Hamiltonian of the spin tube with antiferromagnetic 
legs to a Hamiltonian of spinless fermions using Jordan- Wigner transformation, then map 
it to a modified Luttinger model. The bosonic expressions for spin operators are 



5..$A + — e'(""/") sin ( Va^^x) 
na V / 



where $ and 6 are the bosonic dual fields, and x is defined on the lattice, xj = ja, a is a 
short- distance cutoff. 



(32) 



The bosonized form of the Hamiltonian of the non-interacting chains is 

1 



u 



'Hx = - I dx 



(33) 



where Il\{x) = d^Qx is canonically conjugate momentum to $a- The Luttinger liquid 
parameters are fixed from the Bethe ansatz solution^- 

TtVI - A2 



K 



TT 



(34) 



2(7r — arccos A) 2 arccos A 

The velocity u vanishes and K diverges for A = — 1. This corresponds to the ferromagnetic 
instability point of a single chain. 

The interchain interactions fl30|) between the nearest neighbor chains reads as 



,AA' 



'^'^ (5.$a)(4$a')+^7i [ 7^ COS (v^($, + 



7r 



(27ra)2 
dx 



/dx f \ f* doc 

p-^coB (V4^(*. -*,.))+ 33 y p^cos (v/?(e. - e,.)) 
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+94 
+ 95 



cos (v^(eA - Oy)) COS (^V47r($A + $ 
COS (v^(0A - ey )) COS (^V4^($a " $ 



(35) 



where 5^1 = — 2J| , 5(2 = 2J| , 5(3 = 27rJ^^^y , 94 = 95 = The HamUtonian ( 13T|1 

of the next-to-nearest couphngs has a similar form with a formal change J^. — ?■ J'^.i, 

_L,AA _L,AA 

9i^ 91 = -2J^^^y, 92^ 92 = 2-^l^AA' 

Following the route of Refs.— 1^ it is convenient to introduce a symmetric mode and 
three antisymmetric ones $^25 "^as 



= ^ (<^'i + $2 + $3 + '^'4) , 

"^a, = |($1 + $2 - $3 - '^'4) , 
•^a^ = |($1 - $2 - $3 + '^'4) , 
= i ($1 - $2 + $3 - $4) . 



(36) 



■ as 2 

In terms of the new fields the quadratic part of the Hamiltonian fl29l) is diagonalized to 

3 



u 



"Hn = — i dx 



1 



i=l 



1 



with 



u 



13 



M (^1 + 


2KJI 


un 




= u(^l 


n ( 1 - 


2KJI 


UTT 



+ 



Ks = K 1 + 



UTT 



(37) 



UTT 



■UTT MTT 

KJ\ 



UTT 



+ 



un 



2KJI ^ KJ, 



UTT 



un 



(38) 



The relevant and marginally relevant terms of the interchain couplings are given by 

2 



H 



int 



1=1 



dx 



{2na) 



cos ( V47r$, ) cos f v47r$ 



-291 



dx 



{2nay 



cos v47r<l>, cos v47r<l>„. 



(27ra)2 
dx 



cos V47r$a^ cos \/47r$a3 +2^: 



+2c/3 J (27ra)2 (VTrOag) + 2^3 y 



dx 



{2nay 

+94 



(27ra)2 

dx 
(2^ 



cos ^V^<l>ai j COS ^y^Vr^aa) 

COS (VnQai) COS (v/vrQaa) 



^^^^p [cos (\/7r(ea2 + Bag)) COS 74^(^5 + ^ai^ 
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+ cos (0F(ea, - Gag)) COS (^y4^($s " ^aj) + COS {y^{Qa2 " ©a,,)) COS (^Vi^l^s - ^aj) 

+ COS (v^(eai + Oag)) COS (V^{(l>s + $a: 



(39) 



+^1 / (2^1)2 (v^(0ai + ©aa)) COS (^V^(<l>, + $, 

+ COS (0F(e„, - e„j) COS (y47l{<^, -'^as) 

The terms are irrelevant and are omitted. 

The Hamiltonian (l37Il describes four independent gapless spin- 1/2 chains coupled by the 
interchain interaction in the form of Eg. (1391) . It is expected that the interleg coupling results 
in the Haldane gap in the excitation spectrum. Note that in the vicinity of the single chain 
ferromagnetic instability, A = — 1, the effective bandwidth collapses, m — )■ 0, and the effect 
of the interleg couplings becomes crucial. To find a detail behavior of the gap in the phase 
diagram with antiferromagnetic ( J_l, J'j_ > 0) interleg coupling and ferromagnetic leg regime 
(A < 0) we use the renormalization group analysis. 

The RG equations are derived through the standard technique (see Ref.~, for example). 
The result is 



dg2 
dl 



dg'2 
dl 



[2-2Ka,]g'^, 



dgs 
dl 



1 

2-4 



1 1 

+ 



K K 



93, 



dl 



2 - 



2K„ 



93^ 



dgi 
dl 



dg'i 
dl 



dKs 
dl 



-2(71 



0.3 
2 



2K, 



ai 



9i, 



29I 



Ks 



9i 



Ks 



20 



1/ 93 \ f 93 \ , ^ f 9^ \ ,1/ 94 



dl 2^' \27TUaJ ^^y2TTUaJ 4^^ \27TUa, 



+ (7^1 +Ut^] ■ (40) 



47rMa3/ 2 v47rMa3 

One sees that the gi terms are relevant for Kg + Ka^ < 2; the g[ term is relevant for 
Kg + < 2; the g2 terms are relevant for Kag + < 2; the g'2 term is relevant for 
Ka^ < 1; the g^ term is relevant for K~_^ + < 8; the g'^^ term is relevant for > 1/4. 
Despite the (74 and g'^ terms are irrelevant they are the most relevant terms which couple 
the symmetric and antisymmetric modes^-. 

Using the RG equations the behavior of the gap in the whole phase diagram can be 
established. Following to standard routine, we analyze the effect of the transversal (J^) 
and the longitudinal ( parts of the interleg coupling separately. 

1. Transversal part of the interleg interactions 

In this case J^, J_^^ 7^ and JJ, J^^ = 0, the initial values of the coupling constants are 
given by g,{l = 0) = g[{l = 0) = Q, g^il = 0) = g',{l = 0) = 0, g^il = 0) = 27r J^, g',{l = 0) = 
2ttJ_^^, gi{l = 0) = TxJ^ and g'jj = 0) = txJ^'^. The bare Luttinger parameters are Us = u, 
Ua, = Ua, = Ua, = u, and K,{1 = 0) = K, KaAl = 0) = Ka,{l = 0) = KaS = 0) = K . 
The term g^ is relevant for —1 < A < while the g^ term is irrelevant. It is easily checked 
numerically that g^, g'^ grow whereas g^, g'^ decrease under the RG. It means that 0i, G2, O3 
are locked in one of the vacuum states (9i = 82 = 0, O3 = ^/vr or 61 = ©2 = y/TT, 63 = 
provided >J_^^), fluctuations of the fields 9^. {i = 1,2,3) are completely suppressed. 
Therefore arbitrary > J^^ 7^ generate a gap in the antisymmetric modes (Gj are 
pinned, disordered). 

After the fluctuations of the fields Bj are stopped, the infrared behavior of the symmetric 
mode is governed by the term of the coupling with the antisymmetric modes 

Uint = 2^4 1 cos( ^4^$,) cos(y4^$„,) 
(27ra)^ 

+ 2^4 f — ^cos(V4^$,)cos(V4^<l>a3). (41) 
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where 

94. = MCOS {y/TT[Qa^ + ©ag])) = ^?4(cOS {y/TT[Qa2 + ©ag])), 
94 = ^4(C0S (0r(©ai + ©aa))), 

and (^4, are renormalized couplings provided 93,9'^ = 0{1). Here, the invariance of T^int 
given by Eq.([39D under ©„, ^ -©„, yields (cos (0F(©a, - ©a,))) = (cos (0F(©a, + ©a,)))- 
Despite e**"« has exponentially decaying correlations due to the ©a- are pinned, a scrupulous 
analysis^^ shows that the effective Hamiltonian for $s presents a standard sine-Gordon 
Hamiltonian 

HeS = ^ j dx 

where Kg is a renormalized value of K, and 9 is a. new effective coupling constant. From 
the correlation function (exp [z\/167r$s(x)] exp [iVm^^siy)]) = {ay\x-y\Y it follows 
that the 9 term has a scale dimension AKg. Therefore, it is relevant for Kg < 1/2, when $s 
is pinned, i.e. becomes massiv^. 

To summarize, the transversal part of the interleg coupling supports gapped antisymmet- 
ric modes, the symmetric sector is gapped at Kg < 1/2, and remains gapless at > 1/2. 
The condition Kg = 1/2 determines a boundary between the gapless Spin Liquid XYl 
phas^, and a generalization of the gapped Rung Singlets phase^I for the for-leg spin tube 
(Fig JTOl) . (Hereinafter, we retain names of phases used in the theory of spin ladders with 
ferromagnetic legs.) In the last case, spins on the same rungs, or along the shortest diagonals 
form singlet pairs by a dynamical way. 

2. Longitudinal part of the interleg interactions 

For the case of the longitudinal part of the interleg exchange, J^, J_^^ = and J]_, j'^ 7^ 
0, the bare values of the coupling constants are given by 9i{l = 0) = — 2J|, 5f^(/ = 0) = 
-2Jl^ 92{l = 0) = 2J1, g^{l = 0) = 2J;^ 9^{l = 0) = g'S = 0) = 0, and g,{l = 0) = g'S = 
0) = 0. 

The strong-coupling phase diagram in the vicinity of ferromagnetic instability point (A = 
— 1 and J|, j'^ = 0) obtained in the RG analysis is shown in Fig. [TTl In the sector denoted as 
spin liquid II phas^ the gi^2, 5'i 2 terms are irrelevant. The symmetric and antisymmetric 
modes remain gapless. In the sector marked as a Haldane phase the terms 5^1^25 5^1,2 



Ac 



^9 



dx 



(27ra) 



cos 



(42) 
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FIG. 10: The ground-state phase diagram in the vicinity A = — 1 of the four-leg tube with 
transverse couphng between legs. 

relevant. Since all of the modes are coupled and locked together, both the symmetric and 
antisymmetric modes are gapped. 

The phase of a ferromagnet with antiphase interchain order arises as a result of the 
ferromagnetic instability with increasing interleg antiferromagnetic coupling. The boundary 
of the transition into the phase is obtained by studying the velocity renormalization of the 
corresponding gapless excitations. We mark the transition at Mq- = (i = 1, 2, 3). 

3. Isotropic interleg exchange 

The initial values of the coupling constants are gi(l = 0) = —2J±, g[(l = 0) = —2j'j_, 
g^{l = 0) = 2Jx, 92{l = 0) = 2j;, gS = 0) = 27rJx, g'S = 0) = 27rJl, gS = 0) = vrJ^ and 
g^{l = 0) = ttJ^. 

From the RG equations (HOl) it is seen that the most relevant operators are the g^, g'^ terms. 
Therefore, the antisymmetric sector is gapped, and Bq- are locked in the disordered phase. 
As in the case of the transversal interleg interactions an effective sine-Gordon Hamiltonian 
for the symmetric mode determines phase boundary Kg = 1/2 between gapped and gapless 
phases. Numerical analysis shows that the ground state phase diagram consists of the 
disordered Rung Singlet gapfuU phase and the stripe ferromagnetic phase with dominating 
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FIG. 11: The ground-state phase diagram in the vicinity A = — 1 of the four-leg tube with 
longitudinal coupling between legs. 

intraleg ferromagnetic ordering. The sector of the rung singlet phase increases at — )■ J±_ 
(see Fig. [12]). 

A possible physical picture that reconciles the phase diagram with the results of the 
previous sections might look as follows. The system BIPNNBNO is an array of loosely 
coupled ferrimagnetic chains in a presence of an extremely weak XY anisotropy and a strong 
frustration, ~ J^, what corresponds to the region of the disordered Rung Singlet phase 
close the point of the ferromagnetic instability, A = — 1. 



V. CONCLUSION 



We have studied magnetization process in two-dimensional compound BIPNNBNO which 
exhibits ferrimagnetism of non-Lieb-Mattis type. The investigation is complicated by a lack 
of reliable information about exchange interactions in the system. For a start, we proposed 
the naive model of non-interacting ferrimagnetic chains and showed that an appearance both 
1/3 and 2/3 plateaus can be explained within the model. This provides us the intrachain 
exchange couplings Jaf and Ji. By setting these parameters in the exact diagonalization 
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FIG. 12: The ground-state phase diagram in the vicinity A = — 1 of the four-leg tube with 
isotropic coupling between legs. The phase boundary between the disordered rung singlet gapfuh 
phase and the stripe ferromagnetic phase is shown by the dotted and solid lines for = 0.2 

and J'^/Ji_ = 1.0, respectively. 

routine the magnetization curve for the 32 and 40-sites clusters is numerically calculated. 
We demonstrate that the magnetization curve similar to that observed in the experiment 
is obtained in the regime of weak interchain coupling, Ji ^ J^n > Jnnn- Another revealed 
phenomenon is that a width of the singlet plateau increases with a growth of the antiferro- 
magnetic frustrating coupling between the next-to-nearest chains. Following these results, 
we apply on the tube lattice two low-energy theories which could explain an appearance of 
the singlet phase. The first one is based on the effective XXZ Heisenberg model in a longitu- 
dinal magnetic field in the limit where the interchain coupling dominates, Jnn > ^nnn ^ Ji- 
We derive the critical field destroying the singlet plateau, it turns out that it does not 
depend on the frustration parameter Jnnn- Another analytical strategy is realized via the 
abelian bosonization formalism which is relevant for the opposite limit Ji ^ Jnn > Jnnn- 
We demonstrate that the gapfull disordered Rung Singlets phase comes up when the XY 
exchange anisotropy may tilt the balance from the long-range order with an antiphase in- 
terchain arrangement of ferrimagnetic chains towards the spin liquid phase. A role of the 
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anisotropy in a formation of the spin gap in the original two-dimensional system deserves a 
further study. 
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Appendix A: Wave functions of the ring 



The states of the ring |5'i2<S'345'M) are composed from the functions |mim2m3m4), where 
mj = ±1/2 marks the spin- 1/2 state of the separate i-th (1,1/2) block in the ring 
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The basic functions of the singlet states are given by 
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The triplet states read as 
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Appendix B: Saturation field in the limit of the strong ring coupling 

To get the saturation field the functions of the ring with the total spins S — h and S — ^ 
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3 

with the energies = Jaf/2 — Jnn/3 + Jnnn/2, Eq = 2Jaf + Jnn + 2Jnnn are needed. 

By introducing the pseudo-spin operators in the restricted space similar to Eq. fl23|) the 
original spin operators are presented as follows {S = 1, s = 1/2) 



--) 

22/, 



S. 



S;. 



^in 24 12 



12 ^ 6 



(Bl) 



The effective XXZ spin chain Hamiltonian has the form ( 122|) with the parameters J^y 



-2Ji/3, Jf = Ji/18 and Bf = -7Ji/9 + B - B^, where B^ = Eq- Eg. 

By performing a Jordan- Wigner transformation followed by a particle-hole 
transformation^^ the new Hamiltonian of spinless holes is 



N 



N 



N 



n=l 



n=l 



n=l 



where t = V = Jf and fin = Jf - Bf. 

The saturation field Ssat corresponds to the chemical potential where the hole band starts 
to fill up, fih = -2t. This yields B,^t = B^ + Ji/6 = 3 Jaf/2 + 47^^/3 + 3 Jnnn/2 + Ji/6. 
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